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METHOD OF DISCRETE MODELING AND ITS APPLICATION 
TO ESTIMATION OF TF30 ENGINE VARIABLES 
by Dale J. Arpasi and William M. Bruton 
Lewis Research Center 

SUMMARY 

A method of discrete modeling is presented that effectively isolates steady -state 
model accuracy from dynamic model accuracy. It is formiilated so as to allow the in- 
clusion of readily available steady -state engine relations to achieve any degree of steady- 
state modeling accuracy desired. A dynamic model based on the impulse transfer func- 
tion may then be incorporated without affecting steady-state modeling accuracy. 

The dynamic model is generated by applying a step disturbance of a manipulated var- 
iable to an open-loop engine simvilation. The response of a particular engine variable of 
interest is sampled at a fixed interval to form a set of weighting factors that represent 
the dynamic response of the variable about an operating point. These factors are then 
used to accumulatively weight the past values of the derivative of the manipulated var- 
iable. The summation is multiplied by a simple trim function to account for dynamic 
variations due to changes in operating point. The resultant dynamic model is summed 
with the steady -state model to form the complete discrete model. 

The model was applied to the estimation of the corrected low-pressure-compressor 
rotor speed, turbine-inlet temperature, and corrected airflow of both compressors in the 
TF30-P-3 turbofan engine. These estimates were formed for three different configura- 
tions of the steady -state model by using the derivative of corrected main -combustor fuel 
flow as the argument of the dynamic model. The dynamic model was generated from 
weighting factors obtained at a single operating point. A straight-line function of high- 
pressure -compressor rotor speed was chosen to modify the model dynamics for opera- 
tion at other operating points. Each estimate was evaluated against its actual counter- 
part in the engine simulation for a step in throttle position from a subidle setting to the 
full augmentation setting of power -level position while the simulation was operating under 
the engine's bill -of -material control system. 
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INTRODUCTION 

Airbreathir^ engine control requires the close regulation of certain engine variables 
to provide performance without violating engine operational limits. This regulation can 
be accomplished directly by the controller if these variables can be measured. However, 
the significant variables are not directly measiirable. Sensors for these desired meas- 
urements are often unavailable for either technical or economic considerations. Conse- 
quently, the control engineer is forced to devise either a method of indirect regulation of 
these variables or an accurate means of estimating these variables to permit direct reg- 
ulation in the control. 

Presently, because accurate estimates are difficult to obtain, indirect regulation of 
unmeasurable variables is the common practice. Indirect regulation is normally ac- 
complished by imposing limits on the manipulated variables (e. g. , the acceleration 
limit schedule used in the main fuel controls). These limits are not easily defined in 
terms of the actual limit requirements of the unmeasurable variables (e. g. , compressor 
surge and turbine overtemperature limits) and also are not easily modified for engine-to- 
engine variation and engine deterioration. 

Variable estimation as an alternative to indirect regulation is enhanced by the trend 
toward incorporating sophisticated digital computers into the engine control loop. The 
compvitational power and memory capability of the computer permits storage of a math- 
ematical model of the engine. This model would be capable of providing indirect obser- 
vations of important engine variables based on the actual measurement of other engine 
variables. The indirect observation or estimate of these variables would then be incor- 
porated with some control law for direct regulation. 

This report describes a method of mathematical modeling that is particularly well 
suited to digital engine control. The method is based on the impulse transfer fvinction 
and is formulated completely in the discrete time domain. Consequently, the model may 
be formed empirically by a digital computer directly monitoring an engine simvilation or 
a well-instrumented prototype of the engine. In this case the models are determined by 
a digital computer from an extensive hybrid simulation of the engine. After a discussion 
of the general mathematical formulation of the model and certain practical considera- 
tions, the method of modeling is demonstrated by using it to provide estimates of various 
variables in the TF30-P-3 turbofan engine. 


TF30-P-3 ENGINE DESCRIPTION 

The Pratt & Whitney TF30-P-3 er^ne (fig. 1) is an axial, mixed-flow, augmented, 
twin-spool, low -bypass -ratio turbofan. A single inlet is used for both the fan airflow and 
the engine core airflow. Airflow leaving the fan is separated into two flow streams; one 
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stream passing through the engine core, and the other stream passing through the an- 
nular fan duct. A nine -stage, low-pressure compressor (LPC) includes the core portion 
of the three-stage fan and is connected by a throughshaft to the three-stage, low- 
pressure turbine. A seven-stage, high-pressure compressor (HPC) is connected by a 
hollow shaft to the single-stage, high-pressure turbine. A seventh -stage LPC bleed sys- 
tem discharges air into the fan duct at flight Mach numbers above 1.75. A twelfth -stage 
HPC bleed system automatically discharges air into the fan duct for compressor surge 
control. In the combustor, airflow from the compressors reacts with the injected fuel, 
producing high-energy gas to power the turbines. 

The engine core and fan duct flow streams combine in an augmentor and are dis- 
charged through a variable -area convergent nozzle. The augmentor consists of a dif- 
fuser section and a combustion chamber. The fuel is introduced into the combustion 
chamber through spray rings that are arranged into five separate zones. 

The bill-of -material control (BOM) of the TF30-P-3 engine is used to provide tran- 
sient evaluation of the estimator models and is described in reference 1. The control 
provides regulation of the main -combustor fuel flow, the twelfth-stage compressor bleed, 
the variable exhaust nozzle, and the fuel flow to the five augmentor zones. 


SIMULATION AND ESTIMATION APPARATUS 

The simulation of the TF30-P-3 engine is detailed in reference 2. The equations 
describing the engine were implemented by using the Lewis Research Center’s Analog 
and Hybrid Computing Facility. The Lewis EAI Model 690 Hybrid Computing System was 
used in conjunction with two EAI 231 -R Analog Computers. 

The bill -of -material control, the formulation procedixre of the estimator models, 
and the estimator models were programmed on a general-purpose digital control system 
also located in this facility. This system uses an SEL 810B digital computer and is de- 
scribed in reference 3. Table I summarizes the capabilities of the digital control 
system. 


MATHEMATICAL MODEL FORMULATION 


In applying impulse response to model development, it was assumed that the engine 
is a single -input, multiple -output system with the relation between input and ou^t being 
governed by external conditions. It was also asstimed that measurement of the engine 


variable to be modeled is available (as from a simulation of the ei^ine } even inougn 
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actuality that variable is unmeastmable. 
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Accuracy problems involved with direct application of the impulse transfer function 
to model development are indicated and thus an alternative form, making use of a sep- 
arately evolved steady -state model, is presented. 


Basic Model 

Assume that figure 2(a) represents a source of engine data having single input u(t) 
and containii^ two variables of interest: x(t), which is measurable, and y(t), which is 
unmeasxirable in the actual engine. Let v(t) represent the external conditions of the en- 
gine that modify the relations between the engine variables x(t),y(t) and the manipulated 
variable u(t). Assume that these relations are time -invariant but nonlinear functions of 
u(t) as well as of v(t), but for constant v(t) = Vj^. they become linear for small disturb- 
ances of u(t). 

If u(t) is to be generated by a digital controller, samplers must be incorporated, as 
shown in figure 2(b). The variables u(n), x(n), and y(n) describe u(t), x(t), and y(t), 
respectively, only at the sampling instants. Define the system to be at rest at the n - 1 
sampling instant and u(n - 1) = u^^, x(n - 1) = x^^, and y(n - 1) = y^. If u(n) = u^^ and 
u(n + i) = u^, where i is any positive integer greater than zero, fee system can be de- 
fined in terms of a weighting sequence if the resultant approximate impulse (shown in 
fig. 2(c)) of magnitude Ujj “ ^a small enough so as to maintain the simulation in the 
linear region about x and y„ . The resultant discrete impulse transfer function 
(ref. 4) 

OO 

- ya = E gi("a’Va^[u(n - i) " ^a] 

i=0 

completely describes y(t) around u^^ and v^^ at the sampling instants. The weighting 
factors %(u^)Vg^) are determined from the sampled impulse response and given by the 
determination equation 




y(n + i) - y^ 


b a 


where i = 1 — <» 


( 2 ) 


To evaluate fee steady -state accuracy of the model given by equation (1), let u(n) be 
adjusted so feat u(n - i) = u^ for all i values. The variable y(n) will achieve steady 
state at 
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(3) 


OO 

y(n) > Ya + E 8i(“a' - "a' 

Note that the accuracy of this steady-state model depends upon the accuracy with which 
the values of gj have been determined. 

The accuracy of the weighting factors is primarily a function of three considerations: 

(1) Inherent accuracy of the data source 

(2) Measurement accuracy 

(3) Adherence to the linear region of the model 

If the data source is a simulation, errors will arise as a result of the inability to accu- 
rately define effective areas, volumes, and friction. If the data source is an engine, 
subtle differences in the fabrications of similar engines could cause weighting factor 
error. 

Meastmement accuracy problems are not easily avoided. The impulse disturbance 
of an engine transmits little power to the data source and could easily lead to signal -to- 
noise ratios that would significantly affect the g. accuracy. The primary consideration 
affecting gj accuracy, however, is the necessity of strict adherence to the linear region 
in model applications. Should u be such that the operating -point moves outside this 
region, the g. may become extremely inaccurate and may result in gross errors in both 
the d 3 mamic and steady states. 

Steady -state nonlinearities can be easily handled by generating appropriate functional 
relations amoi^ the various engine variables. These relations are usually available 
from the engine design program and are easily refined on actual engines by usii^ steady- 
state sensors, which are generally available. For these reasons it is not prudent to tie 
steady-state accimacy to dynamically obtained weighting factors. 

To isolate steady-state accuracy from dynamic accuracy, it is necessary that in the 
steady state the model reduce to 


y(n) =fy[u(n),x(n),v(n)] (4) 

where the right side of this equation represents an accurate steady -state model of the en- 
gine. This is easily accomplished by changing the data somrce configuration given in fig- 
ure 2(b) to that given in figure 2(d). In this figure, the dummy variables u' and y' are 
defined as 


u'(n) = u(n) - u(n - 1) 


(5) 


y’(n) =y(n) - fy[u(n),x(n), v(n)] 


( 6 ) 
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If an impulse is applied at u'(n), just as described for the development of the weighting 
factors of equation (2), a set of weighting factors r^ may be determined from 


[•i(u. 




y'(n + i) - Ya 
%-"a 


( 7 ) 


These weighting factors then define the following discrete impulse transfer function: 


y'(n) -y^ = E 

i=0 

By definition, y' and u’ are zero because they represent steady-state values. Asa 
result of substituting equations (5) and (6) into equations (7) and (8) and generalizing to 
any operating point specified by Uj and Vj^, the determination equation and the model 
equation become 


ri(u.,vj^) 


y(n + i) - fy[u(n + i),x(n + i),VjJ 


Ui. - U- 

b 


(9) 


oo 

y(n) = fy[u(n),x(n),v(n)J + E rj(u^, Vj^) [u(n - i) - u(n - i - 1)] (10) 

In steady state, u(n + i) - u(n + i - 1) = 0, and consequently the model of eqtiation (10) 
reduces to that specified by equation (4). This effectively isolates steady -state accuracy 
from the weighting factors and consequently from dynamic acciiracy. Note that in the 
development of the weighting factors of equation (9), the impulse is applied to u’(n). The 
actual disturbance of the manipulated variable is, therefore, a step function. 


Practical Considerations 

The model given in equation (10) is not practical. To build this model, the infinite 
summation must be made finite. This is possible if the engine does not contain a free 
integrator (i. e. , y(n) eventually reaches steady state for constant u(n)). This has al- 
ready been implicitly assmned in the discussion of steady-state relations. Equation (10) 
can then be rewritten as 
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( 11 ) 


y(n) = £y [u(n),x(n),v(n)] + J ^i(^j>\) ['^(” " D - u(n - i - 1)] 

where m is determined by the resolution of the digital controller to be that interval at 
which « 0 for all i > m. 

The amount of information necessary to implement equation (11) depends on the en- 
gine dynamics and nonlinearities and on the desired accuracy of the y(n) estimate. 
Steady-state accuracy is assxired if fy can be defined over the operational envelope of 
the engine (for all v(t)) and can be contained within the storage capability of the control. 
Dynamic accuracy is not only a function of the storage capability of the control but also a 
function of its calculation -time capability. The calculation time of the estimate is essen- 
tially proportional to m. The time may be reduced by extending the control update in- 
terval, thereby reducing m. Extension of the update interval also allows for longer cal- 
culation times but increases the amount of vmdesirable digital dynamics inserted by the 
control into the closed -loop system. 

Since dynamic simulation precision may be questionable, the weighting factor method 
of dynamic identification (or any d 3 mamic identification method based on a simulation) is 
also questionable. It has been indicated that slight inaccuracies in the determination of 
the steady-state relations could also compromise the dynamic accuracy of equation (11). 
Dynamic accuracy may be further compromised because of the dynamic nonlinearity of 
the engine and the inability to completely define the weighting factors for all possible 
values of v(t). Even if adequate definition were possible, storage requirements of the 
weighting factors could be prohibitive . 

To alleviate these problems, the next step in the model formulation is the incorpora- 
tion of a trim function for the model d 3 mamics. Since the steady -state and dynamic 
terms of equation (11) are essentially distinct, the d 3 mamic trim is easily incorporated 
without affecting steady -state accuracy, as follows: 

m 

y(n) = fy [u(n),x(n) ,v(n)J + T^. [x(n),v(n)J £ ri(uj, Vj^) [u(n - i) - u(n - i - l)j (12) 

where T^, is the trim function. The arguments of T^., x(n), and v(n) imply that it could 
vary as a function of the engine operating point. Investigation of this function was made 
empirically and is discussed in a subsequent section of this report. 
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TF30-P-3 ESTIMATOR MODELS 


By using equation (12), models were formulated for estimating the TF30-P-3 engine 
variables given in table n. Turbine-inlet temperature T^ and LPC rotor speed Nj^ 
were corrected to fan-face conditions (station 2 in fig. 1). The airflows and 

Wa jj were corrected to their respective compressor -face conditions (stations 2. 1 and 
2. 2,’ respectively). To form a model using equation (12), three elements are necessary: 

(1) Steady -state model 

(2) Dynamic weighting factors 

(3) Trim function 

In discussing the assumptions and simplifications used in generating each element, no 
attempt is made to optimize the model accuracy by var 3 dng the element complexity. A 
very basic approach is taken in the generation of each element in order to obtain a worst - 
case evaluation of the modeling technique. 


Steady -State Model Generation 


It is apparent from equation (9) that the weighting factors and consequently the dy- 
namic accuracy of the model depend heavily on the argument of the steady -state model. 
The argument that will provide the most accurate dynamic model will be that which min- 
imizes the absolute magnitude of the weighting factors in the dynamic term of the model. 
This follows from the assumption that the dynamic accuracy of the model is more ques- 
tionable then the steady-state accuracy. For this reason, three steady-state models 
were considered, each with arguments of significantly different dynamic content. These 
arguments (CW^, CPg g 2’ "^g^ given in table m and were corrected to the fan- 

face conditions (station 2 in fig. 1) . Gage temperature T is defined as 

o 


T 


g 


1 

1 + 0* 05s 


(T3 - T2 + Tg) 


(13) 


Manufacturer’s data indicate a good dynamic correlation between Tg and T^. 

The steady-state functions fy were derived from the hybrid engine simtilation by 
techniques described in reference 5. A slow ramp in CW^ was inserted into the sim- 
vilation while the exhaixst nozzle was held at minimum area (its normal position during 
subaugmentation operation). The minimum and maximum CWj settings were deter- 
mined by the simulation's range of operation at the two flight conditions necessary for 
choked and unchoked nozzle operation. The flight conditions used were 

(1) Sea-level static (unchoked nozzle operation) 

(2) Mach number, M, 1.2; altitude, H, 12 192 meters (choked nozzle operation) 
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The minimum and maximum CWj were 0. 113 and 0,997 kilogram per second, respec- 
tively. At each flight condition a ramp of 100-second duration was used to ensure 
steady -state engine operation. 

The sampled values of the measurable and unmeasurable engine variables were 
stored in tables of 55 words each at 0. 0227-kilogram-per-second increments of CWj. 
Each table required nine words of heading. Tables for each variable were formed for 
three bleed conditions (closed, twelfth-stage open, and seventh-stage open) and at the 
two flight conditions. The storage requirement for the steady-state functional relation 
of each variable was 3 x 2 x 64 = 384 words. 

The resultant function tables provide the steady -state relations of the engine varia- 
bles (whether measurable, x(n), or unmeasurable, y(n)) to the manipulated variable 
CWj(n) at the sample instant n. Note that, if the engine is not in steady state, these 
relations will provide the values that the engine variables would eventually achieve if 
CWj were held at CWj(n), Assume that f^ and fy^ represent the functional relation 
between x(n) and CWj(n) and between y(n) and CWj(n), respectively, and that f^^ is a 
monotonically increasing function. Then the inverse function f”^ unique and f , the 
steady-state fimctional relation between y(n) and x(n), may be determined by using the 
equation 


The steady -state functional relations between each variable and CWj are given in 
figure 3. Only closed compressor bleed conditions are given for the sake of simplicity. 
The effects of compressor bleed are discussed in reference 5. Note that the curves rep- 
resenting unchoked nozzle operation blend with those representing choked nozzle opera- 
tion at high values of CWj. This is due to the nozzle choking at high rotor speeds even 
at sea-level-static conditions. Also note that the final discrete model does not depend on 
this method of steady-state modelii^. Any desirable method of steady -state modeling 
may be used. 


Weighting Factor Generation 

The determination equation (eq. (9)) was used to generate the weighting factors for 
the mathematical models. The rj(u.,Vj^) were determined for the Vj^ given in table IV 
for each Uj given in table V. When the simulation achieved steady state at CW^ = Uj, 
CWj was stepped to a value of tij + 0. 113 kilogram per second. This magnitude was 
chosen as the minimum to give good measurement resolution. Each y-variable and x- 
variable was sampled, corrected as indicated previously, and stored. The sample time 
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was chosen to be 50 milliseconds to allow adequate calculation time for the estimates. 
According to reference 1, this is the maximum interval that may be used without dynamic 
degradation. Sixty -four samples were taken of each variable per CWj disturbance at 
each condition, allowing for a total transient time of 3. 2 seconds. 

The sampled data contained sufficient information when coupled with the steady-state 
functions of figure 3 to solve equation (9) for the weighting factors. Weighting factors 
could then be developed by using either CWj or the transient data for x(n) as the argu- 
ment of the steady-state function relating it to y(n). The result at each sample instant 
was s\ibtracted from the sampled y(n) at that instant. The r.(u.,Vjj) were formed at 
each sample i by dividing the result by the step magnitude. 

From equation (9), specific weighting factors for the CWj input were developed as 




^ y(n + i) - fyw[cWj(n + i)' 


0.113 


(15) 


at Vj^ = Vj^ and for the u. in table V. These are given in figure 4 to illustrate the non- 
linearities of the weighting factors as a function of setpoint Uj. The range of vij repre- 
sents a range in CNjj of 9788 to 14 230 rpm, which is greater than the range in CNjj 
from idle to intermediate. 

Since the argument of the steady-state function used to obtain these weighting factors 
is CWj, the weighting factors represent the response of each variable to a step in CW^ 
at each operating point. The i = l factor is approximately the gain of this response. 
This is true since in the first 50 milliseconds the measured y(n + 1) has barely moved, 
whereas the steady-state model has reached its final value. The time for the response 
to reach zero is proportional to the time constant; and overshoot, if any, is proportional 
to system damping. Note that, as j increases, the time to reach steady state de- 
creases. The difference between the final value of the weighting factors and zero (if 
steady state has been reached) represents inaccuracies in the steady-state model. 

It should be pointed out that the step in CWj occurred at i = 0. Since the engine is 
nonanticipator y, the weighting factor at i = 0 is zero. 

The effect of HPC bleed on the weighting factors is also illustrated in figure 4. At 
Uj = Uj^, the HPC bleed was opened and the u^ = u^^ weightii^ factors were again gen- 
erated. As compared with the effects of variations of Uj, the bleed effects on the 
weighting factors are negligible. 

Figure 5 compares the effects of setpoint u^ and LPC bleed at the external condi- 
tion Vjj = Vg on weighting factors again generated by using equation (15). The effect of 
LPC bleed is also negligible. The effects of setpoint, however, are notably different 
from those at Vj^ = Vj (fig. 4). The difference is due to exhaust nozzle choking effects. 
At Vj^ = Vj (fig. 4) the exhaust nozzle is unchoked for setpoints u^,Ug, and Ug. At 
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corresponding setpoints at the e:diaust nozzle is choked. The major difference 

between the weighting factors of figure 4 and those of figure 5 occurs at these setpoints. 
F^vme 6 shows weighting factors developed from the determination equation 


■•yP.th’V ' 


y(n + 1 ) - 2. 2^” * '*]} 


0.113 


(16) 


at Vj^ = and for the Uj in table V. Here the inverse functional relation of equa- 
tion (14) has been used to relate the unmeasurable variables directly to CPg 2 2 * 
in this case the argvunent of the steady -state model is CPg 2 2 

equation (15), these weighting factors reflect the relative d^amics of the y-variables to 
CPs 2 2 closer they are dynamically, the smaller the absolute magnitude of the 
we^hting factors). Note that CNj^ and CW^ are dynamically very close to 

^^s,2.2- 

In f^ure 7, weighting factors were generated by using equation (17). As expected, 
CT^ is dynamically very close to CT^: 


■^yT,l<“j>V = 


y(n + i) - f, 


* *'1l 


0.113 


Storage and timing requirements for the dynamic model depend upon the number of 
sets of weighting factors needed for each variable and also upon the number of weighting 
factors contained in each set. For the weightii^ factors shown in figures 4 and 5, 64 
words of data storage are necessary per condition per variable. The computation of the 
dynamic model requires a minimum of 64 multiplications and additions and 128 loads and 
stores. For the digital computer described in table I this amounts to a minimum of 
0. 561 millisecond of computation time. 


Trim Function Generation 

Review of figures 4 and 5 makes it apparent that the problems associated with de- 
veloping weighting factors to exactly define engine dynamics for all possible steps in 
CW| at all ei^ine operating points are prohibitive. It would be desirable to be able to 
incorporate a simple trim function with a single set of weighting factors that would allow 
adequate definition of engine dynamics vinder all circxmi stances. Also, since the dy- 
namics of the computerized engine simulation from which the weighting factors wovild be 
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developed are assumed to contain some degree of inaccviracy, a trim function may be 
necessary in applying the estimator models to actual engine control systems. 

Instead of a detailed analysis to determine an optimum trim function for each varia- 
ble, an experimental generation of a simple trim function was performed. The trim 
function obtained by using this method provided a worst -case evaluation of the discrete 
estimator model. A more detailed analysis in the generation of the trim function could 
only improve the accuracy of the estimator model. 

The first step in the experimental development of a trim function was to choose a 
single set of weightii^ factors for use in the estimator model. The weighting factors 
correspondii^ to choked nozzle operation were chosen since under the majority of cir- 
cumstances the engine is operating under this condition. The weighting factors obtained 
at lowest CNjj are seen from figures 5 to 7 to contain the greatest dynamic content. 

The weighting factors for use in the estimator models, therefore, were chosen to be 
those represented by the solid lines of figures 5 to 7. 

The second step in the development was to choose the argument of the trim function. 
The purpose of the trim function is to modify the weighting factors to account for the var- 
iation in engine dynamics from that encountered at the point at which the weighting factors 
were obtained. A large contribution to the dynamic variation can be related to the mech- 
anical speed of the rotors (fig. 4). For this reason the dynamic trim was made a func- 
tion of mechanical speed. For simplicity, a straight-line function was chosen, the inter- 
cept and slope of which were empirically derived. 

The trim function Tj,(Njj) was developed experimentally for the estimation of 
CW^ at Vj^ = Vg . The estimator model for CW^ was chosen as 


CW 


A, L " VlW 


(CWf) 


64 P 

Tj.(Nh) Z - ^Wf(" - ' - 


(18) 


This was done by observing the error between the measured value of CW^ ^ 
estimated value for a 4-second ramp in CWj from 0, 113 kilogram per second to 0. 793 
kilogram per second. A straight-line function was chosen, and its intercept and slope 
were varied until a minimum error response for the ramp was observed. The resultant 
function was 


Tj.(Njj) = -1.25X10"^ Njj + 1. 875 (19) 

The use of Njj in the model dictates that it must be measurable. 

The solid line in figure 8 is a plot of the trim function. The optimum trim function 
should have some relation to the area contained in the weighting factor curves of figure 5, 
relative to the weighting factors obtained at Uj^ and V2 for each variable. As an 
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attempt to verify this, a relative dynamic power content (RDPC) was formed as 


RDPCj = 



^yW,i^"j’'^2^ 


1/2 


( 20 ) 


The RDPC's for selected variables are shown in figure 8 for comparison with the empir- 
ically derived trim function. 

Note that the RDPC's are closer to the trim function at h^h values of Njj. This 
seems to indicate that matchir^ of dynamics is more significant at high rotor speeds than 
at low speeds. Note also that the RDPC's of CW^ fall very close to the line, suggest- 
ing a definite correlation between the relative dynamic power and the optimum trim func- 
tion. No attempt was made to further optimize the trim function or to find a better cor- 
relation with the optimum. Further work in this area could only improve the results 
presented herein, and these are considered sufficient to evaluate the concept of discrete 
modeling. 


EVALUATION OF ESTIMATOR MODELS 

Three estimator models were evaluated in this study. These were 

64 

y(n) = fyw(CWf) + T^(Njj) E ryw,i("pV 2 )[cWj(n - i) - CWj(n - i - 1)] (21) 

64 

y(n) = fyw[^Pw(C^s,2.2^] + 2 ^yP,iS»’^ 2 ^E^Wj(n - i) - CWj(n - i - 1)] (22) 

64 

y(n) = fyw[f-yCTg)] + T^(Njj) E ryT^i(^i,V 2 )[cWf(n - i) - CWj(n - i - 1)] (23) 

where the y(n) are given in table II, the ry^ ^ are given in figure 5 (u^ = u^, bleed 
closed), the ryp j are given in figure 6 (uj = Uj), and the ry,j, j are given in figure 7 
(Uj = Uj^). The steady -state functions are given in figure 3. 

The estimates developed from equation (21) for the variables given in table n by us- 
ing the trim fvmction described in equation (19) are compared with those estimates ob- 
tained by using a \mity trim function T^(Njj) = 1 in figure 9. The transients shown in 


13 


figure 9 and in figures 11 to 13 represent the difference between the actual measurement 
of the variable on the hybrid simulation and its estimated value over a typical engine 
transient. The results are shown at the external operating coiulitions of table IV. 

The engine transient used in generating these data was a throttle step from idle to 
full axigmentation with the engine simulation operating under its bill-of -material control 
(ref. 1). The transients used to obtain the data of figures 9 and 11 to 13 are given in fig- 
ure 10 in terms of fuel flow rate Wj for each of the Vj^ in table IV. 

The value of the trim function is apparent in figure 9 by comparing the solid line 
(trim defined by eq. (19)) to the dashed line (unity trim). At nozzle choked conditions 
(V 2 ,V 3 , and v^), estimation errors are generally smaller than those occurrii^ at an vin- 
choked nozzle condition (v ^). This may be due in part to the trim function being devel- 
oped at Vg but is due mostly to the weightii^ factors used in the model for dynamic es- 
timation at condition Vj also being obtained at V£. (Note the difference between the 
weighting factors of fig. 4 and those of fig. 5. ) 

Improvement in the accuracy of the estimates of figvme 9 might be expected with the 
use of a more complex dynamic model. However, improvement could also be obtained 
with a model of equivalent complexity if the argument of the steady -state model were 
made more dynamically accurate. The more dynamically accurate this argmnent be- 
comes, the smaller the weighting and consequently the smaller the effects of operating 
condition on variable dynamics. 

Figure 11 compares the estimates developed from equation (22), which uses 
CPs 2 2^® argument of the steady-state term, for v^^ to v^. And figure 12 com- 
pares the estimates developed from equation (23), which uses CTg as the argument of 
the steady-state term, for to v^. Figure 11 shows definite improvement in estima- 
tion accuracy over figure 9 for all variables, with the exception of CT^ and possibly 
CWa jj. This indicates a closer dynamic correspondence of CW^ L 
CPs 2 2 between these variables and CW^. 

Figure 12 indicates a strong dynamic correlation between CT^ and CTg. This fol- 
lows from the large reduction in CT^ estimation error over that resulting from the 
models of equations (21) and (22). Note that the final error in each transient of this fig- 
ure is significantly larger than the final error in the transients shown in figures 9 and 11. 
This implies inaccuracies in the steady -state model relatii^ CTg to CWj. 

Airflow and combustor temperature are perhaps the two most unmeasurable varia- 
bles in any engine. The preceding discussion suggests that CPg 2 2 ^"^g 

possibly be used without dynamic compensation as substitutes for CW^ jj and CT^, 
respectively. To evaluate this possibility, the following estimator models were formed: 

^^A, H ^ ^WHw[^PW^^^s, 2. 2^] 
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(25) 


CT 


4 



which relate high compressor corrected airflow and corrected combustor temperature to 
CPs 2 2 ^®spectively, through the appropriate steady -state relations. The 

results are presented in figure 13. Although the estimation errors are somewhat larger 
than those obtained with dynamic compensation, these estimates could still prove useful 
in some applications. 


CONCLUSIONS 

The method of discrete modeling described in this report provides not only a 
straightforward means of generating estimates of unobservable variables, but also pro- 
vides the control engineer with a useful means of evaluating relations between engine 
variables. The method is easily automated. It is tailored for formulation and imple- 
mentation on a digital computer, allowing for easy adaptation into a digital control 
system. 

Isolating the steady -state model from the dynamic model permits the original 
steady-state design equations, or any simplification of these equations, to be included 
into the model without consideration of whether the dynamics to be included affect steady- 
state accuracy. Dynamic model accuracy on the other hand may be enhanced though the 
proper choice of the argument of the steady -state function. The closer dynamically the 
observed argument is to the variable to be estimated, the smaller the dependence of the 
model on the weightii^ factors and consequently the higher the model accuracy. 

Applying the modeling method to the estimation of TF30-P-3 engine variables re- 
veals that compressor bleeds have little effect upon dynamic accuracy as compared with 
the effects of engine setpoint and degree of exhaust nozzle choking. The effects of engine 
setpoint chaises may be somewhat accounted for by applyii^ a straight-line function of 
high-pressure -compressor (HPC) rotor speed as a trim of a dynamic mode developed at 
a single setpoint. 

Finally, it may be concluded that, with a proper steady -state relation, gage temper- 
ature (as described herein) and low -pres sure -compressor (LPC) discharge static pres- 
sure may be used as alternate measurements for turbine-inlet temperature and HPC cor- 
rected airflow, respectively. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 14, 1976, 

505-05. 
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APPENDIX - SYMBOLS 


CN corrected rotor speed, rpm 

2 

CP corrected pressure, kg/m 

CT corrected temperature , K 

CW corrected flow rate, kg/sec 

CWj corrected fuel flow rate, kg/ sec 

f functional relation 

g generalized we^hting factor 

H altitude , m 

i sample indexii^ variable (integer) 

M Mach number 

N rotor speed, rpm 

n integer 

2 

P pressTxre, kg/m 

RDPC relative dynamic power content 
r dynamic weighting factor 

s Laplace operator, sec“^ 

T temperature, K 

Tj. trim function 

t time, sec 

u generalized manipulated variable 

Uj generalized internal engine operating condition 

Vj^ generalized external engine operating condition 

W flow rate, kg/sec 

Wj fuel flow rate, kg/sec 

X generalized observable variable 

y generalized unobservable variable 

z z -transform operator 
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Subscripts: 

A 

a,b,c 

e 

g 

H 

i 

L 

m 

PW 

s 

ss 

TW 

T4W 

WHW 

WLW 

WNW 

XU 

yP 

yT 

yu 

yw 

yx 

1 , 2 , 2 . 1 , 2 . 2 , 
3,4,4. 1,5,6 

Superscript: 


air 

particular values 

estimate 

gage 

high-pressure compressor 
sample indexing variable 
low-pressure compressor 
measured 

functional designation relating g 2 
static 

steady state 

functional designation relating CT^ to CW^ 
functional designation relating CT^ to CWj 
functional designation relating CW^ jj to CWj 
functional designation relating CW^ to CWj 
functional designation relating CNj^ to CW^ 
functional designation relating x to u 
functional designation relating y to CP 
functional designation relating y to CT 
functional designation relating y to u 
functional designation relating y to CWj 
functional designation relating y to x 

ei^ine station numbers as defined in fig. 1 
difference variable 
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TABLE I. - DIGITAL SYSTEM CAPABILITIES 


Digital computer 

Magnetic core memory size, words 

16 384 

Word length, bits plus parity 

16 

Memory cycle time, nsec 

750 

Add time, (isec 

1.5 

Subtract time, jusec 

1.5 

Multiply time, jusec 

4.5 

Divide time, jusec 

8.25 

Load time, jitsec 

1.5 

Store time, /isec 

1.5 

Indirect addressing 

Infinite 

Indexing 

Total memory 

Priority interrupts 

28 Separate levels 

Index registers: 


Independent 

1 

In conjunction with lower accumulator 

1 

Physical size, cm (in.); 


Width 

60. 1 (24) 

Height 

157. 5 (62) 

Depth 

76. 2 (30) 

Interval timers 

Complement 

2 

Accuracy, clock pxilses 

±1 

Clock rates, kHz 

572 , 286, 160, 143, 80, 71.5, 


40, 35.75, 20, 10 

Counter 

16 -Bit binary 

Output 

Priority interrupt to computer 

Analog acquisition unit 

Number of multiplexers, digitizers, samples, and holds 

2 

Overall sample rate (maximum), kHz 

40 

Resolution of digital data, bits 

12 (plus sign) 

Output code 

Two’s complement 

Number of channels 

64 

Input range, V full scale 

±10 

Input impedance, MJ2 (shunted by 10 pF) 

10 

Maximum source resistance, 

1000 

Conversion time, jusec 

38 

Input setting time, jusec 

9 

Sample -and-hold aperture time, nsec 

500 

Safe input voltages, V 

±20 sustained 


±100 for less than 100 jusec 

Total error with calibration, percent 

0.073 



TABLE I. - Concluded. DIGITAL SYSTEM CAPABILITIES 


Frequency acquisition unit 

Nximber of channels 

10 

Nature of input 

Continuously varying or pulsatile 

Resolution of digital data, bits 

12 

Switch selectable clock rates, kHz 

20, 80, 100, 400, external 

Overall accuracy, bits 

±1 

Undate rate 

Once per cycle of input frequency 

Maximum input frequency, kHz 

1 

Input amplitude range 

100 mV to 30 V peak to peak 

Analog output unit 

Total number of digital -to -analog 

26 

conversion channels 


Resolution (10 channels), bits 

12 (plus sign) 

Resolution (16 channels), bits 

11 (plus sign) 

Output voltage range, V full scale 

±10 

Output current (maximum), mA 

10 

Output impedance, 

<1 

Accuracy (12 bit), percent of full scale 

±0. 1 

Accuracy (13 bit), percent of full scale 

±0.05 

Slew rate, V/jusec 

1 

Settling time for 10-V step to within 0. 05 

20 

percent of final value, /isec 


Logical output unit 

Number of electronic switch outputs 

32 

Number of contact closure outputs 

32 

Maximum voltage, V 

30 

Maximum current, mA 

100 

Priority interrupt processor 

Number of channels 

10 

Input impedance, kJ2 

47 

Input voltage range, V 

±10 

Comparator switching 

Trigger on rise or fall 

Comparator hysteresis 

Adjustable from 35 mV to 650 mV 

Comparator output, V 

+7 

Monostable multivibrator: 


Pulse width, jisec 

3 

Pulse height, V 

+7 



TABLE n. - DESCRIPTION OF ENGINE VARIABLES 


TO BE ESTIMATED 


Estimated 

variable 

Description 

CNl 

Corrected low-pressure -compressor rotor 
speed, Nl/|/^ 

CWa,L 

Corrected low-pressure -compressor airflow, 


^A,2.2^®2.1^/ }/h. 1 

CWa,h 

Corrected high-pressure -compressor airflow, 


CT4 

Corrected turbine-inlet temperature, T^/02 


TABLE m. - DESCRIPTION OF ARGUMENTS TO BE 


USED IN THE GENERATION OF ESTIMATES 


Argument 

variable 

Description 

CWj 

^^s,2.2 

CTg 

Corrected main-combustor fuel flow, 

Corrected low -pressure -compressor discharge 
static pressure, Pg g 2^^ 

Corrected gage temperature, Tg/^2 


TABLE IV. - EXTERNAL ENGINE TEST CONDITIONS 


Generalized external 
engine operating 
condition, 

Vj^, where k = 

Engine external conditions 

Mach number, 
M 

Altitude, 

H, 

m 

1 

0 

0 (sea-level static) 

2 

1.2 

12 192 

3 

1.2 

3 048 

4 

2.2 

15 240 





TABLE V. - ENGINE SETPOINTS FOR 


WEIGHTING FACTOR 
DETERMINATION 


Generalized internal 
engine operating 
conditions, 

Uj, where j = 

Corrected main -combustor 
fuel flow (CW^) setpoint, 
kg/ sec 

1 

0. 113 

2 

.227 

3 

.340 

4 

.453 

5 i 

, 680 


/-Fan inside 
/ diameter 
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(b) Discretized version of the simplified block diagram. 



Sample, i 


(c) Approximate impulse with duration of 1 sample period. 

u(n) 



(d) Dynamic model development. 


Figure 2. * Mathematical model development. 
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Corrected LPC discharge static pressure, Corrected LPC speed, CNl, rpm Corrected turbine-injet temperature. 

CPc ? 9, kg/m^ 



(e) Corrected low-pressure -com pressor (LPC) discharge static- Corrected gage temperature function, 

pressure function, fpy^. 


Figure 3. - Steady-state relations for choked and unchoked nozzle conditions. 






(38s;6)|)/(D9S/6)1) 

■(lA‘!n)WHMj ‘jopej Su!)i|6j3M moijjib odH papsjjoo 
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Sample, i 

(b) Corrected high -pressure -compressor (HPC) airflow [y(n) ^ ^(n)]. 

Figure 4. - Effects of cor rected-fuel -flow-rate CW| operating point and high -pressure-compressor bleed on weighting factors derived by using CWj as the argument of the 
steady-state model and for a 0. 113-kilogram -per-second step in CWj at the generalized external operating condition V|^ = vj. 












Generalized internal 
engine operating 
condition 

j » u^. bleeds closed 
j = U2, bleeds closed 
j = U3, bleeds closed 
j = U4, bleeds closed 
j = U5, bleeds closed 


(D8s/ 6)|)/)| ‘jopej 

6 uuMb|dM djn)ejaduia) laiuj-eujqjn) papajjoo 


(Das/6)|}/(Das/6 )|) ‘(^a 
‘jope; buuLjbjaw moip!b odH papaJJOQ 
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Sample, i 

(b) Corrected high-pressure-compressor (HPC) airflow [y(n) = j^(n)]. 

Figure 5 - Effects of corrected-fuel-flow-rate CWf operating point and low-pressure-compressor bleed on weighting factors derived by using corrected gage 
temperature CTq as the argument of the steady-state model and for a 0 . 113 -ki log ram -per -second step in CW^ at the generalized external operating condi- 


Generalized interna! 
engine operating 
condition 

ij ' Uj, bleeds closed 
ij = U2. bleeds closed 
ij = U3, bleeds closed 
I; = u>i. bleeds closed 



V ^ (D0S/6)|)/ (305/6)0 ‘(^A‘fn)^^^j 

(30S/6)|)/ujdj *(2A‘fn)MNj 'paads Qj-j pepajjoQ ‘jope^ 6uuM6|aMMO|iJiB Odl papajJOD 
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Sample, i 

(d) Corrected low -pressure-compressor (LPC) airflow [y(n) = CW^ |_(n)]. 
Figure 5 . - Concluded. 







Generalized internal 
engine operating 
condition 



(d) Corrected low-pressure-compressor (LPC) airflow [y(n) = L(n)]. 

Figure 6. - Concluded. 




Generalized internal 
engine operating 
condition 
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Sample, i 

(b) Corrected high -pressure-compressor (HPC) airflow [y(n) * CW^ 

Figure 7. - Effects of corrected-fuel-flow-rate CWj operating point on weighting factors derived by using corrected gage temperature CTg as the argument of 
the steady-state model and for a 0. 113-ki log ram -per -second step In CWj at the generalized external operating condition with both bleeds closed. 












Corrected fuel flow, 
CWf, kg/ sec 


1.5r- 



Time, sec 


(a) = vj. (b) V|^ = V2- (c) v,^ = V3. (d) v,^ = V4. 

Figure 10 , - Fuel-flow-rate Wf transients due to a throttle step from idle to full augmentation with 
engine operating under bill-of-material control at selected external conditions V|^. 






Corrected LPC airflow Corrected turbine-inlet 

estimation error. temperature estimation 

(CWa, l)^ - (C^A, Og. - (CT4)g. 

kq/sec K 



gage temperature CTg. 



Time, sec 


(b-1) V[( * Vj. (b-2) Vj( = V 2 - (b-3) V|^ = V 3 . (b-4) = v^. 

(b) Corrected low-pressure-compressor (LPC) airflow estimation error as function of corrected low- 
pressure-compressor (LPC) discharge static pressure CP^^ 2 . 2- 

Figure 13. - Estimation errors for steady-state models at selected external conditions V|^. 
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